clear all
clc
x=sym('x')
f1=sqrt(x+x^3)-(x+x^2)^(1/3)
taylor(f1,x,4)
